Alleviation of arsenic-induced neurobehavioral defects with selenium in the larvae of Zaprionus indianus.

Selenium is an essential antioxidative micronutrient. This study was conducted to characterize the arsenic toxicity induced on the African fig fly, Zaprionus indianus, and its possible amelioration by selenium. We used computational tools and in vivo experiments to elucidate the mechanism of action of arsenic and selenium on Z. indianus larvae. We conducted experiments to study neurobehavioral parameters including learning and memory ability test and crawling and contraction assays. Our in silico study revealed twelve primary targets of arsenic trioxide. The gene ontology annotation of primary and secondary targets of arsenic trioxide revealed selenocysteine metabolic processes as one of the most reliable targets. To validate our in silico data, we analyzed the effect of arsenic trioxide on larvae of Z. indianus and tested the possible amelioration by sodium selenite supplementation. Our data demonstrated that the arsenic trioxide deteriorated the learning and memory ability of 2nd instar larvae of Z. indianus and such effect was reversed by sodium selenite supplementation. Furthermore, crawling and contraction assay done on 3rd instar larvae showed that there was reduction in both parameters upon arsenic trioxide exposure, which was restored with sodium selenite supplementation. Altogether, our computational and in vivo results strongly indicated that the neurobehavioral defects induced by arsenic trioxide on the larvae of Z. indianus can be successfully alleviated in the presence of sodium selenite.

A computational study on mitogenome-encoded proteins of Pavo cristatus and Pavo muticus identifies key genetic variations with functional implications.

BACKGROUND: The Pavo cristatus population, native to the Indian subcontinent, is thriving well in India. However, the Pavo muticus population, native to the tropical forests of Southeast Asia, has reduced drastically and has been categorised as an endangered group. To understand the probable genetic factors associated with the decline of P. muticus, we compared the mitogenome-encoded proteins (13 proteins) between these two species. RESULTS: Our data revealed that the most frequent variant between these two species was mtND1, which had an alteration in 9.57% residues, followed by mtND5 and mtATP6. We extended our study on the rest of the proteins and observed that cytochrome c oxidase subunits 1, 2, and 3 do not have any change. The 3-dimensional structure of all 13 proteins was modeled using the Phyre2 programme. Our data show that most of the proteins are alpha helical, and the variations observed in P. muticus reside on the surface of the respective proteins. The effect of variation on protein function was also predicted, and our results show that amino acid substitution in mtND1 at 14 sites could be deleterious. Similarly, destabilising changes were observed in mtND1, 2, 3, 4, 5, and 6 and mtATP6-8 due to amino acid substitution in P. muticus. Furthermore, protein disorder scores were considerably altered in mtND1, 2, and 5 of P. muticus. CONCLUSIONS: The results presented here strongly suggest that variations in mitogenome-encoded proteins of P. cristatus and P. muticus may alter their structure and functions. Subsequently, these variations could alter energy production and may correlate with the decline in the population of P. muticus.

An immunoinformatics approach to study the epitopes of SARS-CoV-2 helicase, Nsp13.

Introduction and objective: Vaccines are administered worldwide to control on-going coronavirus disease-19 (COVID-19) pandemic caused by SARS-CoV-2. Vaccine efficacy is largely contributed by the epitopes present on the viral proteins and their alteration might help emerging variants to escape host immune surveillance. Therefore, this study was designed to study SARS-CoV-2 Nsp13 protein, its epitopes and evolution. Methods: Clustal Omega was used to identify mutations in Nsp13 protein. Secondary structure and disorder score was predicted by CFSSP and PONDR-VSL2 webservers. Protein stability was predicted by DynaMut webserver. B cell epitopes were predicted by IEDB DiscoTope 2.0 tools and their 3D structures were represented by discovery studio. Antigenicity and allergenicity of epitopes were predicted by Vaxijen2.0 and AllergenFPv.1.0. Physiochemical properties of epitopes were predicted by Toxinpred, HLP webserver tool. Results: Our data revealed 182 mutations in Nsp13 among Indian SARS-CoV-2 isolates, which were characterised by secondary structure and per-residue disorderness, stability and dynamicity predictions. To correlate the functional impact of these mutations, we characterised the most prominent B cell and T cell epitopes contributed by Nsp13. Our data revealed twenty-one epitopes, which exhibited antigenicity, stability and interactions with MHC class-I and class-II molecules. Subsequently, the physiochemical properties of these epitopes were analysed. Furthermore, eighteen mutations reside in these Nsp13 epitopes. Conclusions: We report appearance of eighteen mutations in the predicted twenty-one epitopes of Nsp13. Among these, at least seven epitopes closely matches with the functionally validated epitopes. Altogether, our study shows the pattern of evolution of Nsp13 epitopes and their probable implications.

Introducción y objetivo: Las vacunas se administran a nivel mundial para controlar la pandemia en curso de la enfermedad por coronavirus de 2019 (COVID-19) causada por SARS-CoV-2. A la eficacia de la vacuna contribuyen ampliamente los epítopes presentes en las proteínas virales, y su alteración puede contribuir a que las variantes emergentes se escapen de la vigilancia inmunológica del huésped. Por tanto, este estudio fue diseñado para estudiar la proteína Nsp13 de SARS-CoV-2, sus epítopes y su evolución. Métodos: Se utilizó Clustal Omega para identificar las mutaciones de la proteína Nsp13. La estructura secundaria y la tasa de desorden se predijeron mediante los servidores web CFSSP y PONDR-VSL2. La estabilidad de la proteína fue predicha mediante el servidor web DynaMut. Los epítopes de las células B fueron predichos mediante las herramientas DiscoTope 2.0 de IEDB, y sus estructuras en 3D fueron representadas mediante Discovery Studio.La antigenicidad y alergenicidad de los epítopes fueron predichas mediante Vaxijen2.0 y AlergenFPv.1.0. Las propiedades fisioquímicas de los epítopes fueron predichas mediante Toxinpred, la herramienta del servidor web HLP. Resultados: Nuestros datos revelaron 182 mutaciones en Nsp13 entre los aislados indios de SARS-CoV-2, que fueron caracterizadas mediante las predicciones de la estructura secundaria y la capacidad de desorden por residuo, la estabilidad y la dinamicidad. Para correlacionar el impacto funcional de estas mutaciones, caracterizamos los epítopes más prominentes de las células B y las células T a los que contribuyó Nsp13. Nuestros datos revelaron veintiún epítopes, que exhibieron antigenicidad, estabilidad e interacciones con las moléculas MHC de clase I y clase II. Seguidamente se analizaron las propiedades fisioquímicas de estos epítopes. Además, en estos epítopes de Nsp13 residen ocho mutaciones. Conclusiones: Reportamos el aspecto de ocho mutaciones en los veintiún epítopes de Nsp13 predichos. Entre estos, al menos siete epítopes concuerdan estrechamente con los epítopes funcionalmente validados. En su conjunto, nuestro estudio refleja el patrón evolutivo de los epítopes de Nsp13 y sus implicaciones probables.

RNA degradation eliminates developmental transcripts during murine embryonic stem cell differentiation via CAPRIN1-XRN2.

Embryonic stem cells (ESCs) are self-renewing and pluripotent. In recent years, factors that control pluripotency, mostly nuclear, have been identified. To identify non-nuclear regulators of ESCs, we screened an endogenously labeled fluorescent fusion-protein library in mouse ESCs. One of the more compelling hits was the cell-cycle-associated protein 1 (CAPRIN1). CAPRIN1 knockout had little effect in ESCs, but it significantly altered differentiation and gene expression programs. Using RIP-seq and SLAM-seq, we found that CAPRIN1 associates with, and promotes the degradation of, thousands of RNA transcripts. CAPRIN1 interactome identified XRN2 as the likely ribonuclease. Upon early ESC differentiation, XRN2 is located in the nucleus and colocalizes with CAPRIN1 in small RNA granules in a CAPRIN1-dependent manner. We propose that CAPRIN1 regulates an RNA degradation pathway operating during early ESC differentiation, thus eliminating undesired spuriously transcribed transcripts in ESCs.

Single-cell RNA-seq reveals early heterogeneity during aging in yeast.

The budding yeast Saccharomyces cerevisiae (S. cerevisiae) has relatively short lifespan and is genetically tractable, making it a widely used model organism in aging research. Here, we carried out a systematic and quantitative investigation of yeast aging with single-cell resolution through transcriptomic sequencing. We optimized a single-cell RNA sequencing (scRNA-seq) protocol to quantitatively study the whole transcriptome profiles of single yeast cells at different ages, finding increased cell-to-cell transcriptional variability during aging. The single-cell transcriptome analysis also highlighted key biological processes or cellular components, including oxidation-reduction process, oxidative stress response (OSR), translation, ribosome biogenesis and mitochondrion that underlie aging in yeast. We uncovered a molecular marker of FIT3, indicating the early heterogeneity during aging in yeast. We also analyzed the regulation of transcription factors and further characterized the distinctive temporal regulation of the OSR by YAP1 and proteasome activity by RPN4 during aging in yeast. Overall, our data profoundly reveal early heterogeneity during aging in yeast and shed light on the aging dynamics at the single cell level.

Non-uniform aspects of the SARS-CoV-2 intraspecies evolution reopen question of its origin.

Several hypotheses have been presented on the origin of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) from its identification as the agent causing the current coronavirus disease 19 (COVID-19) pandemic. So far, no solid evidence has been found to support any hypothesis on the origin of this virus, and the issue continue to resurface over and over again. Here we have unfolded a pattern of distribution of several mutations in the SARS-CoV-2 proteins in 24 geo-locations across different continents. The results showed an evenly uneven distribution of the unique protein variants, distinct mutations, unique frequency of common conserved residues, and mutational residues across these 24 geo-locations. Furthermore, ample mutations were identified in the evolutionarily conserved invariant regions in the SARS-CoV-2 proteins across almost all geo-locations studied. This pattern of mutations potentially breaches the law of evolutionary conserved functional units of the beta-coronavirus genus. These mutations may lead to several novel SARS-CoV-2 variants with a high degree of transmissibility and virulence. A thorough investigation on the origin and characteristics of SARS-CoV-2 needs to be conducted in the interest of science and for the preparation of meeting the challenges of potential future pandemics.

An issue of concern: unique truncated ORF8 protein variants of SARS-CoV-2.

Open reading frame 8 (ORF8) shows one of the highest levels of variability among accessory proteins in Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of Coronavirus Disease 2019 (COVID-19). It was previously reported that the ORF8 protein inhibits the presentation of viral antigens by the major histocompatibility complex class I (MHC-I), which interacts with host factors involved in pulmonary inflammation. The ORF8 protein assists SARS-CoV-2 in evading immunity and plays a role in SARS-CoV-2 replication. Among many contributing mutations, Q27STOP, a mutation in the ORF8 protein, defines the B.1.1.7 lineage of SARS-CoV-2, engendering the second wave of COVID-19. In the present study, 47 unique truncated ORF8 proteins (T-ORF8) with the Q27STOP mutations were identified among 49,055 available B.1.1.7 SARS-CoV-2 sequences. The results show that only one of the 47 T-ORF8 variants spread to over 57 geo-locations in North America, and other continents, which include Africa, Asia, Europe and South America. Based on various quantitative features, such as amino acid homology, polar/non-polar sequence homology, Shannon entropy conservation, and other physicochemical properties of all specific 47 T-ORF8 protein variants, nine possible T-ORF8 unique variants were defined. The question as to whether T-ORF8 variants function similarly to the wild type ORF8 is yet to be investigated. A positive response to the question could exacerbate future COVID-19 waves, necessitating severe containment measures.

The importance of accessory protein variants in the pathogenicity of SARS-CoV-2.

The coronavirus disease 2019 (COVID-19) is caused by the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS- CoV-2) with an estimated fatality rate of less than 1%. The SARS-CoV-2 accessory proteins ORF3a, ORF6, ORF7a, ORF7b, ORF8, and ORF10 possess putative functions to manipulate host immune mechanisms. These involve interferons, which appear as a consensus function, immune signaling receptor NLRP3 (NLR family pyrin domain-containing 3) inflammasome, and inflammatory cytokines such as interleukin 1ß (IL-1ß) and are critical in COVID-19 pathology. Outspread variations of each of the six accessory proteins were observed across six continents of all complete SARS-CoV-2 proteomes based on the data reported before November 2020. A decreasing order of percentage of unique variations in the accessory proteins was determined as ORF3a > ORF8 > ORF7a > ORF6 > ORF10 > ORF7b across all continents. The highest and lowest unique variations of ORF3a were observed in South America and Oceania, respectively. These findings suggest that the wide variations in accessory proteins seem to affect the pathogenicity of SARS-CoV-2.

Periodically aperiodic pattern of SARS-CoV-2 mutations underpins the uncertainty of its origin and evolution.

Various lineages of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) have contributed to prolongation of the Coronavirus Disease 2019 (COVID-19) pandemic. Several non-synonymous mutations in SARS-CoV-2 proteins have generated multiple SARS-CoV-2 variants. In our previous report, we have shown that an evenly uneven distribution of unique protein variants of SARS-CoV-2 is geo-location or demography-specific. However, the correlation between the demographic transmutability of the SARS-CoV-2 infection and mutations in various proteins remains unknown due to hidden symmetry/asymmetry in the occurrence of mutations. This study tracked how these mutations are emerging in SARS-CoV-2 proteins in six model countries and globally. In a geo-location, considering the mutations having a frequency of detection of at least 500 in each SARS-CoV-2 protein, we studied the country-wise percentage of invariant residues. Our data revealed that since October 2020, highly frequent mutations in SARS-CoV-2 have been observed mostly in the Open Reading Frame (ORF) 7b and ORF8, worldwide. No such highly frequent mutations in any of the SARS-CoV-2 proteins were found in the UK, India, and Brazil, which does not correlate with the degree of transmissibility of the virus in India and Brazil. However, we have found a signature that SARS-CoV-2 proteins were evolving at a higher rate, and considering global data, mutations are detected in the majority of the available amino acid locations. Fractal analysis of each protein's normalized factor time series showed a periodically aperiodic emergence of dominant variants for SARS-CoV-2 protein mutations across different countries. It was noticed that certain high-frequency variants have emerged in the last couple of months, and thus the emerging SARS-CoV-2 strains are expected to contain prevalent mutations in the ORF3a, membrane, and ORF8 proteins. In contrast to other beta-coronaviruses, SARS-CoV-2 variants have rapidly emerged based on demographically dependent mutations. Characterization of the periodically aperiodic nature of the demographic spread of SARS-CoV-2 variants in various countries can contribute to the identification of the origin of SARS-CoV-2.

Emerging genetic diversity of SARS-CoV-2 RNA dependent RNA polymerase (RdRp) alters its B-cell epitopes.

The RNA dependent RNA polymerase (RdRp) plays crucial role in virus life cycle by replicating the viral genome. The SARS-CoV-2 is an RNA virus that rapidly spread worldwide and acquired mutations. This study was carried out to identify mutations in RdRp as the SARS-CoV-2 spread in India. We compared 50217 RdRp sequences reported from India with the first reported RdRp sequence from Wuhan, China to identify 223 mutations acquired among Indian isolates. Our protein modelling study revealed that several mutants can potentially alter stability and flexibility of RdRp. We predicted the potential B cell epitopes contributed by RdRp and identified thirty-six linear continuous and twenty-five discontinuous epitopes. Among 223 RdRp mutants, 44% of them localises in the B cell epitopes region. Altogether, this study highlights the need to identify and characterize the variations in RdRp to understand the impact of these mutations on SARS-CoV-2.

Emergence of unique SARS-CoV-2 ORF10 variants and their impact on protein structure and function.

The devastating impact of the ongoing coronavirus disease 2019 (COVID-19) on public health, caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has made targeting the COVID-19 pandemic a top priority in medical research and pharmaceutical development. Surveillance of SARS-CoV-2 mutations is essential for the comprehension of SARS-CoV-2 variant diversity and their impact on virulence and pathogenicity. The SARS-CoV-2 open reading frame 10 (ORF10) protein interacts with multiple human proteins CUL2, ELOB, ELOC, MAP7D1, PPT1, RBX1, THTPA, TIMM8B, and ZYG11B expressed in lung tissue. Mutations and co-occurring mutations in the emerging SARS-CoV-2 ORF10 variants are expected to impact the severity of the virus and its associated consequences. In this article, we highlight 128 single mutations and 35 co-occurring mutations in the unique SARS-CoV-2 ORF10 variants. The possible predicted effects of these mutations and co-occurring mutations on the secondary structure of ORF10 variants and host protein interactomes are presented. The findings highlight the possible effects of mutations and co-occurring mutations on the emerging 140 ORF10 unique variants from secondary structure and intrinsic protein disorder perspectives.

Autoimmunity roots of the thrombotic events after COVID-19 vaccination.

Although vaccination represents the most promising way to stop or contain the coronavirus disease 2019 (COVID-19) pandemic and safety and effectiveness of available vaccines were proven, a small number of individuals who received anti-SARS-CoV-2 vaccines developed a prothrombotic syndrome. Vaccine-induced immune thrombotic thrombocytopenia (VITT) can be triggered by the adenoviral vector-based vaccine, whereas lipid nanoparticle-mRNA-based vaccines can induce rare cases of deep vein thrombosis (DVT). Although the main pathogenic mechanisms behind this rare phenomenon have not yet been identified, both host and vaccine factors might be involved, with pathology at least in part being related to the vaccine-triggered autoimmune reaction. In this review, we are considering some aspects related to pathogenesis, major risk factors, as well as peculiarities of diagnosis and treatment of this rare condition.

Implications derived from S-protein variants of SARS-CoV-2 from six continents.

The spike (S) protein is a critical determinant of the infectivity and antigenicity of SARS-CoV-2. Several mutations in the S protein of SARS-CoV-2 have already been detected, and their effect in immune system evasion and enhanced transmission as a cause of increased morbidity and mortality are being investigated. From pathogenic and epidemiological perspectives, S proteins are of prime interest to researchers. This study focused on the unique variants of S proteins from six continents: Asia, Africa, Europe, Oceania, South America, and North America. In comparison to the other five continents, Africa had the highest percentage of unique S proteins (29.1%). The phylogenetic relationship implies that unique S proteins from North America are significantly different from those of the other five continents. They are most likely to spread to the other geographic locations through international travel or naturally by emerging mutations. It is suggested that restriction of international travel should be considered, and massive vaccination as an utmost measure to combat the spread of the COVID-19 pandemic. It is also further suggested that the efficacy of existing vaccines and future vaccine development must be reviewed with careful scrutiny, and if needed, further re-engineered based on requirements dictated by new emerging S protein variants.

COVID-19 Vaccines and Thrombosis-Roadblock or Dead-End Street?

Two adenovirus-based vaccines, ChAdOx1 nCoV-19 and Ad26.COV2.S, and two mRNA-based vaccines, BNT162b2 and mRNA.1273, have been approved by the European Medicines Agency (EMA), and are invaluable in preventing and reducing the incidence of coronavirus disease-2019 (COVID-19). Recent reports have pointed to thrombosis with associated thrombocytopenia as an adverse effect occurring at a low frequency in some individuals after vaccination. The causes of such events may be related to SARS-CoV-2 spike protein interactions with different C-type lectin receptors, heparan sulfate proteoglycans (HSPGs) and the CD147 receptor, or to different soluble splice variants of the spike protein, adenovirus vector interactions with the CD46 receptor or platelet factor 4 antibodies. Similar findings have been reported for several viral diseases after vaccine administration. In addition, immunological mechanisms elicited by viral vectors related to cellular delivery could play a relevant role in individuals with certain genetic backgrounds. Although rare, the potential COVID-19 vaccine-induced immune thrombotic thrombocytopenia (VITT) requires immediate validation, especially in risk groups, such as the elderly, chronic smokers, and individuals with pre-existing incidences of thrombocytopenia; and if necessary, a reformulation of existing vaccines.

A unique view of SARS-CoV-2 through the lens of ORF8 protein.

Immune evasion is one of the unique characteristics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) attributed to its ORF8 protein. This protein modulates the adaptive host immunity through down-regulation of MHC-1 (Major Histocompatibility Complex) molecules and innate immune responses by surpassing the host's interferon-mediated antiviral response. To understand the host's immune perspective in reference to the ORF8 protein, a comprehensive study of the ORF8 protein and mutations possessed by it have been performed. Chemical and structural properties of ORF8 proteins from different hosts, such as human, bat, and pangolin, suggest that the ORF8 of SARS-CoV-2 is much closer to ORF8 of Bat RaTG13-CoV than to that of Pangolin-CoV. Eighty-seven mutations across unique variants of ORF8 in SARS-CoV-2 can be grouped into four classes based on their predicted effects (Hussain et al., 2021) [1]. Based on the geo-locations and timescale of sample collection, a possible flow of mutations was built. Furthermore, conclusive flows of amalgamation of mutations were found upon sequence similarity analyses and consideration of the amino acid conservation phylogenies. Therefore, this study seeks to highlight the uniqueness of the rapidly evolving SARS-CoV-2 through the ORF8.

Notable sequence homology of the ORF10 protein introspects the architecture of SARS-CoV-2.

The current Coronavirus Disease 19 (COVID-19) pandemic, caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) shows similar pathology to MERS and SARS-CoV, with a current estimated fatality rate of 1.4%. Open reading frame 10 (ORF10) is a unique SARS-CoV-2 accessory protein, which contains eleven cytotoxic T lymphocyte (CTL) epitopes each of nine amino acids in length. Twenty-two unique SARS-CoV-2 ORF10 variants have been identified based on missense mutations found in sequence databases. Some of these mutations are predicted to decrease the stability of ORF10 in silico physicochemical and structural comparative analyses were carried out on SARS-CoV-2 and Pangolin-CoV ORF10 proteins, which share 97.37% amino acid (aa) homology. Though there is a high degree of ORF10 protein similarity of SARS-CoV-2 and Pangolin-CoV, there are differences of these two ORF10 proteins related to their sub-structure (loop/coil region), solubility, antigenicity and shift from strand to coil at aa position 26 (tyrosine). SARS-CoV-2 ORF10, which is apparently expressed in vivo since reactive T cell clones are found in convalescent patients should be monitored for changes which could correlate with the pathogenesis of COVID-19.

Molecular assessment of proteins encoded by the mitochondrial genome of Clarias batrachus and Clarias gariepinus.

The population of catfish, Clarias batrachus has substantially diminished in various countries and studies show that another related species Clarias gariepinus is replacing it. The better adaptability and survivability of C. gariepinus over C. batrachus could be attributed to the metabolic differences between these two species, which is primarily regulated by mitochondrial activities. To understand the reasons behind this phenomenon, we performed in silico analyses to decipher the differences between the proteins encoded by the mitochondrial genome of these two related species. Our analysis revealed that out of thirteen, twelve proteins encoded by the mitochondrial genome of these two species have substantial variations between them. We characterised these variations by analysing their effect on secondary structure, intrinsic disorder predisposition, and functional impact on protein and stability parameters. Our data show that most of the parameters are changing between these two closely related species. Altogether, we demonstrate the molecular insights into the mitochondrial genome-encoded proteins of these two species and predict their effect on protein function and stability that might be helping C. gariepinus to gain survivability better than the C. batrachus.

Urgent Need for Field Surveys of Coronaviruses in Southeast Asia to Understand the SARS-CoV-2 Phylogeny and Risk Assessment for Future Outbreaks.

Phylogenetic analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is focused on a single isolate of bat coronaviruses (bat CoVs) which does not adequately represent genetically related coronaviruses (CoVs) [...].

The molecular assessment of SARS-CoV-2 Nucleocapsid Phosphoprotein variants among Indian isolates.

Coronavirus disease- 2019 (COVID-19) has rapidly become a major threat to humans due to its high infection rate and deaths caused worldwide. This disease is caused by an RNA virus, Severe Acquired Respiratory Syndrome -Corona Virus-2 (SARS-CoV-2). This class of viruses have a high rate of mutation than DNA viruses that enables them to adapt and also evade host immune system. Here, we compared the first known Nucleocapsid Phosphoprotein (N protein) sequence of SARS-CoV-2 from China with the sequences from Indian COVID-19 patients to understand, if this virus is also mutating, as it is spreading to new locations. Our data revealed twenty mutations present among Indian isolates. Out of these, mutation at six positions led to changes in the secondary structure of N protein. Further, we also show that these mutations are primarily destabilising the protein structure. The candidate mutations identified in this study may help to speed up the understanding of variations occurring in SARS-CoV-2.